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Abstract 

The  aerodynamic  section  characteristics  of 
several  helicopter  profiles  have  been  computed 
using  the  Reynolds-averaged,  thin-layer  Navler- 
Stokes  equations,  and  an  eddy-viscosity  model  to 
approximate  boundary-layer  turbulence.  The  calcu¬ 
lations  were  performed  over  a  wide  range  of  Mach 
numbers,  angles  of  attack,  and  Reynolds  numbers, 
with  particular  emphasis  on  transonic  separated 
flow  conditions  that  are  relevant  to  advancing 
olades  in  high-speed  forward  flight.  An  extensive 
study  was  made  of  the  accuracy  of  the  numerical 
code,  ARC2D,  and  of  the  sensitivity  of  the  results 
to  the  numerical  parameters  and  approximations. 
Also,  comparisons  were  made  with  a  large  body  of 
experimental  data.  The  numerical  results  repro¬ 
duce  the  measured  airfoil  behavior  across  the 
transonic  regime,  and  the  details  of  the  computed 
flow  field  provide  new  insights  into  transonic 
airfoil  behavior. 


Symbols 

C,)  «  drag  coefficient 

CQ  ■  transonic  similarity  drag  parameter,  Eq.  (5/ 

«  lift  coefficient 

CM  •  pitching  moment  c. ff lclent  about  x/c  •  0.25 

Cp  ■  pressure  coefficient 

C*  •  sonic  value  of  C_ 

P  P 

c  ■  chord 

E  •  flux  vector,  Eq.  (1b) 
e  •  total  internal  energy  per  unit  volume 
F  •  flux  vector,  Eq.  (1b) 

J  •  Jacobian  of  the  coordinate  transformation 
Ma  •  free-stream  Mach  number 
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M  •  transonic  similarity  Mach  number  parameter, 
Eq.  (6) 

An  •  mesh  spacing  normal  to  the  airfoil  in  physi¬ 
cal  space 

Pr  •  Prandti  number 
p  •  pressure 

Q  -  dependent  variable  vector,  Eq.  (lb) 

Re  ■  Reynolds  number 

S  •  viscous  stress  vector,  Eq.  (1c) 

As  ■  mesh  spacing  along  the  airfoil  In  physical 
space 

t  •  time  In  Eqs.  (1-3);  also  airfoil  thickness  in 
Eqs.  (5,6) 

•  free-stream  velocity 

U  -  contravarlant  velocity  in  the  ^-direction 
u  »  x-component  of  velocity  In  physical  space 
V  -  contravarlant  velocity  in  the  n-dlrectlon 
v  -  y-component  of  velocity  In  physical  space 
x  -  horizontal  coordinate  In  physical  space 
y  •  vertical  coordinate  In  physical  space 
a  ■  angle  of  attack,  Jsg 
T  ■  circulation 
y  •  specific  heat  ratio 

n  •  transformed  coordinate  noimal  tv  the  a <  r 
surface 

k  •  thermal  conductivity 
u  •  viscosity 

Vj  •  turbulent  eddy  viscosity 

i  ■  transformed  coordinate  along  the  airfoil 
surface 
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p  *  density 


Table  i  limitations  of  wind  tunnel  testing 


t  •  time  in  computational  space 


Primary 


Introduction 


The  transonic  flow  phenomena  that  occur  on 
the  advancing  blade  tips  of  modern  helicopters 
have  major  effects  on  the  aerodynamic  performance, 
vibratory  airloads,  and  acoustic  radiation  of  the 
rotor.  It  is  well  known  that  under  transonic 
conditions,  the  aerodynamic  section  characteris¬ 
tics,  such  as  lift,  drag,  and  pitching-moment 
coefficients,  differ  substantially  from  even  the 
qualitative  behavior  of  subsonic  airfoils.  Typi¬ 
cal  nonlinear  behavior  is  illustrated  in  Fig.  1 
from  Ref.  i.  Especially  noteworthy  are  the  vari¬ 
ous  changes  in  the  derivatives  dC^/da  and 
dC„/da  as  a  function  of  Mach  number. 


Wall  interference 
Three-dimensional  effects 
Costs 

Model  construction 
Tunnel  operations 


Secondary 


Free-stream  turbulence 

Nonuni fora  flow 

Reynolds  number  limits 

Geometry  errors  and  model  deformation 

Measurement  errors 

Limited  amount  of  data 

Turn-around  time 


Fluid  Dynamics  (CFD)  technology,  can  be  used  to 
derive  useful,  quantitative  information  about 
transonic  airfoil  characteristics.  While  CFD  is 
free  of  most  of  the  limitations  of  wind  tunnel 
testing,  it  has  its  own  shortcomings  (cf. 

Table  2).  As  a  result  of  the  numerical  approx¬ 
imations,  computations  from  different  numerical 
analyses  are  never  in  perfect  accord,  Just  as 
results  from  different  facilities  seldom  agree. 


Fig.  1  Transonic  Mach  number  characteristics  of 
an  NACA  0012  airfoil.  (From  Ref.  1.) 


Table  2  Limitations  of  computational 
fluid  dynamics 


Primary 


Turbulenoe  model 
Computer  speed  and  memory 
Costs 

Code  and  software 
Computer  operations 


In  general,  airfoil  C^,  CD,  and  Cg  depend 
primarily  upon  the  geometry,  Mach  number,  and 
angle  of  attack;  and  secondarily,  upon  Reynolds 
number,  surface  conditions,  and  free-stream  turbu¬ 
lence.  Heretofore,  a  complete  set  of  nonlinear 
airfoil  characteristics,  such  as  those  shown  in 
Fig.  1 ,  could  only  be  obtained  from  expensive  wind 
tunnel  tests.  However,  the  validity  and  accuracy 
of  airfoil  data  can  be  tainted  by  wall  interfer¬ 
ence.  three-dimensional  effects,  and  other  fac¬ 
tors,  as  indicated  in  Table  1.  As  a  result  of  the 
various  difficulties  and  uncertainties  of 
transonic  airfoil  testing,  results  for  the  same 
model  often  differ  more  from  one  tunnel  to  another 
than  from  one  airfoil  section  to  another.2'4 

This  paper  examines  the  extent  to  which  an 
alternative  tool,  namely,  modern  Computational 


Secondary 


Transition  model 
Computational  grids 
Numerical  errors 
Turn-around  time 


The  approach  in  this  investigation  is  to 
solve  the  time-averaged,  thin-layer  Navler-Stokes 
equations  wltn  a  simple  eddy-viscosity  model  for 
the  turbulent  viscous  layer  next  to  the  airfoil, 
using  the  NASA  Ames  Research  Center  CFD  code 
ARC2D.5'7  To  our  knowledge,  this  is  the  first 
attempt  to  apply  such  a  sophisticated  code  to  a 
wide  range  of  practical  airfoil  cases.  Therefore, 
our  primary  goals  are:  (i)  to  demonstrate  that  we 
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N  \ 


can  reproduce  the  complex  nonlinear  behavior  shoun 
in  Fig.  i,  (2)  to  determine  the  physical  processes 
of  shock-wave  development  and  vtscous-inviscid 
interaction  responsible  for  this  behavior,  and 
!3)  to  reduce  the  uncertainty  surrounding  quanti¬ 
tative  transonic  airfoil  characteristics. 

The  specific  objectives  of  this  investigation 
were  to  obtain  the  lift,  drag,  and  pitching-moment 
behavior  for  a  conventional  airfoil  (e.g.,  the 
NACA  0012)  and  for  a  representative  modern  tip 
section  (the  Vertol  VR-8)  as  functions  of  Mach 
number  for  the  low-lift  conditions  that  correspond 
to  the  advancing  blade  of  a  high-speed  helicop¬ 
ter.  Figure  2  shows  the  combinations  of  Ma  and 
<i  for  which  calculations  have  been  performed. 
Particular  attention  was  given  to  the  transonic 
regime  above  the  drag-divergence  Mach  number; 
namely,  0.8  <  Ma  <  1.1,  where  wind  tunnel  data  are 
more  limited  and  uncertain.  The  results  in  this 
regime  were  supplemented  by  calculations  for 
NACA  0006  and  0015  profiles  to  study  some  effects 
of  airfoil  thickness.  A  limited  number  of  refer¬ 
ence  calculations  were  also  performed  with  the 
invlscid  mode  of  ARC2D  and  with  potential-flow 
codes. 

As  indicated  in  Fig.  2,  we  have  also  calcu¬ 
lated  the  section  characteristics  for  some  typical 


nover  conditions,  i.e.,  0.5  <  <  0.6,  with  par¬ 

ticular  emphasis  on  lift-to-drag  ratios.  On  the 
other  hand,  we  have  avoided  the  conditions  of 
retreating-blade  stall,  because  we  are  skeptical 
that  the  current  state  of  the  art  in  turbulence 
modeling  would  permit  such  cases  to  be  calculated 
with  confidence. 


Description  of  the  Numerical  Method 

The  computer  code  ARC2D  (Refs.  5-7)  developed 
at  NASA  Ames  Research  Center  solves  either  the 
Euler  equations  or  the  thin-layer,  Reynolds- 
averaged  form  of  the  Navler-Stokes  equations  with 
a  relatively  simple  algebraic  turbulence  model, 
for  two-dimensional  steady  or  unsteady  flow.  The 
highlights  of  the  method  are  outlined  below;  the 
details  are  documented  in  Refs.  5  and  6. 

Governing  Equations  and  Boundary  Conditions 

The  two-dimensional  equations  of  motion  for  a 
viscous  perfect  gas  are  written  in  strong  conser¬ 
vation-law  form  in  a  generalized,  body-fitted 
curvilinear  coordinate  system  (C,n,r).  with  the 
thln-shear-layer  assumption  for  high-Reynolds- 
number  turbulent  flows,  as  described  below. 
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and  U  and  V  are  the  contravariant  velocities 
along  the  £  and  n  directions,  given  by 
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The  generalized  coordinate  system 
e  •  e(x.y.t) 
n  -  n(x,y,t) 

T  •  t 


(3) 


allows  the  boundary  surfaces  In  the  physical  plane 
to  be  mapped  Into  rectangular  surfaces  In  the 
computational  domain,  as  Illustrated  In  the  upper 
part  of  Fig.  3.  This  transformation  simplifies 
the  application  of  the  boundary  conditions  In  the 
computational  domain  and  allows  grid  points  to  be 
clustered  in  the  physical  plane  In  regions  of  high 


gradients.  The  metrics  E  ,  etc..  In  Eq.  (1) 

P .  x 


are  formed  from  the  derivatives  of  xT,  x,,  etc., 
and  the  Jacobian  of  the  transformation  given  by 
Eq.  (3). 


The  pressure,  density,  and  velocity  com¬ 
ponents  are  related  to  the  energy  per  unit  volume 
by  the  equation  of  state  for  a  perfect  gas, 
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The  effective  dynamic  viscosity  j  -  j.  *  is 
typically  comprised  of  a  constant  (molecular) 


value, 
ity, 


u^,  plus  a  computed  turbulent  eddy  viscos- 


\ij/ o,  as  discussed  below. 


Zero  velocity  Is  imposed  as  a  boundary  condi¬ 
tion  on  the  body;  the  thermodynamic  variables  at 
the  body  are  derived  from  adiabatic  wall  condi¬ 
tions  and  the  conditions  3p/3n  •  0  and 
3p/3n  *0  at  the  wall.  The  condition  of  uniform 
flow  at  infinity  Is  approximated  at  the  outer 
boundary  of  the  computational  domain  by  a  nonre- 
fleetive  characterlstlc-like  boundary  procedure, 
with  the  additional  stipulation  that  the  veloci¬ 
ties  at  these  outer  boundaries  Include  the  Induc¬ 
tion  effects  of  the  circulation  (or  lift)  of  the 
body.  That  Is,  the  velocity  field  due  to  a  com¬ 
pressible  potential  vortex  of  strength 
r  •  (l/JJU^cC, ,  where  CL  Is  the  lift  generated 
by  the  airfoil,  Is  superimposed  at  the  outer 
boundary  as  a  perturbation  to  the  uniform  free 
stream.  The  effectiveness  of  this  circulation 
correction  In  the  ARC2D  code  was  demonstrated  in 
Ref.  7  for  lnvlscid  calculations,  and  we  verified 
the  process  for  viscous  flows  by  numerical  experi¬ 
ments  In  which  the  distance  to  the  outer  boundary 
was  varied  from  6  to  *8  chords,  within  these 
bounds ,  the  lift  and  drag  were  found  to  be  con¬ 
stant  to  within  0.51  and  3(»  respectively. 


(2)  Approximations  in  the  Viscous  Regions 


At  the  Reynolds  numbers  of  Interest  In  most 
helicopter  applications,  the  effects  of  viscosity 
are  concentrated  near  rigid  boundaries  and  In  wake 
regions.  Furthermore,* the  flow  gradients  In  the 
viscous  layers  of  Interest  are  generally  much 
larger  In  the  direction  normal  to  the  body  than 
they  are  along  the  body.  Therefore,  it  often 
suffices  to  neglect  the  "streamwise"  viscous  terms 
In  the  governing  equations  (see  Eq.  (D).  This 
"thin-layer”  approximation  Is  similar  to  the  clas¬ 
sical  boundary-layer  approximation,  but  the  normal 
momentum  equation  is  retained  and  solved.  Thus, 
the  pressure  Is  allowed  to  vary  across  the  viscous 
layer.  Blottner 1 0  has  recently  discussed  the 
thin-layer  approximation  In  detail. 


The  thin-layer  approximation  can  break  down 
for  low  Reynolds  numbers  and  In  regions  of  massive 
flow  separation.  However,  It  affords  significant 
computational  economy.  The  resultant  equations 
are  slightly  simpler;  but  more  importantly,  grid 
systems  can  be  used  that  concentrate  grid  lines 
near  the  body,  with  fine  grid  spacing  in  the 
normal  direction  and  relatively  coarse  grid  spac¬ 
ing  along  the  body  (see  Fig.  3).  For  the  flow 
conditions  considered  In  this  paper,  we  feel  it  is 
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Fig.  3  Computational  grids. 


far  more  cost  effective  to  neglect  the  streaowise 
viscous  terms  and  to  concentrate  the  grid  along 
the  body  as  shown  in  Fig.  3,  than  it  would  be  to 
solve  the  full  Reynolds-averaged  Navier-Stokes 
equations  with  more  nearly  equal  clustering  of 
grid  points  in  both  the  normal  and  streaowise 
directions. 

The  algebraic  mixing-length  model  of  Baldwin 
and  Lomax11  is  used  to  approximate  the  effects  of 


turbulence  via  the  turbulent  eddy  viscosity,  v^. 

As  discussed  in  Refs,  n  and  12,  the  model  con¬ 
sists  of  an  Inner  layer  that  Is  governed  by  the 
Prandtl  mixing  length  with  van  Driest  damping,  and 
an  outer  layer  that  follows  the  Clauser  approxima¬ 
tion.  Computed  vortlclty  is  used  m  defining  the 
reference  mixing  length.  The  Baldwln-Lomax  model 
was  designed  specifically  for  use  with  the  thin- 
layer  approximation,  and  It  is  likewise  suspect 
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p 


for  massively  separated  flous.  Therefore,  we  have 

not  attempted  to  calculate  C,  ,  for  example. 

max 

The  ARC2D  code  can  be  applied  to  either 
laminar  or  turbulent  flows,  but  transition  to 
turbulence  Is  not  modeled.  Therefore,  transition 
point  locations  on  the  body  have  to  be  specified, 
or  else  determined  by  separate  empiricism,  as  was 
done  In  the  recent  study  by  Mehta,  Chang,  and 
Cebecl12  using  ARC20.  Transition  was  fixed 
aroitrarlly  at  x/c  -  0.01  on  both  the  upper  and 
lower  surfaces  of  the  airfoil  for  all  the  cases 
presented  In  this  paper.  A  limited  study  of  tran¬ 
sonic  cases  with  other  transition  locations  showed 
the  lift  and  pitching  moment  to  be  relatively 
Insensitive  to  this  parameter,  but  the  drag  gener¬ 
ally  decreased  with  Increasing  xT,  as  would  be 
expected. 

Computational  Grids 

Body-conforming  grids  of  both  the  "0"  and  "C" 
topology  were  used  and  evaluated,  with  the 
C-topology  Illustrated  In  Fig.  3  generally  con¬ 
verging  to  grid- Independent  values  with  fewer  grid 
points  or  with  fewer  Iterations.  For  most  of  the 
calculations  shown  In  this  paper,  a  variation  of 
the  hyperbolic  grid-generation  scheme  of  Barth, 
Pulliam,  and  Buning1’  was  used,  typically  with 
193  points  In  the  e  direction  (72  points  on  each 
surface  of  the  airfoil  and  25  In  the  wake)  and  6k 
points  In  the  n  direction.  The  outer  boundary 
was  normally  placed  6  chords  from  the  airfoil. 

The  grid  shown  In  Fig.  3  is  for  illustrative  pur¬ 
poses  only;  the  grid  spacing  shown  near  the  body 
is  much  coarser  than  the  ones  we  actually  used. 

The  location  of  the  outer  boundary  was  found 
to  be  relatively  unimportant  when  the  circulation 
correction  was  used,  as  noted  earlier.  However, 
the  grid  spacing  near  the  body  Is  Important, 
especially  with  regard  to  the  drag  coefficient. 
Whereas  we  found  the  lift  to  be  essentially  Inde¬ 
pendent  of  the  normal  spacing,  An,  for  minimum 
values  less  than  about  0.000k,  we  were  obliged  to 
reduce  Analn  to  about  0.00002  for  reasonably 
reliable  drag  results. 

The  dependence  on  the  spacing  along  the  body, 
As,  was  less  pronounced,  but  it  was  necessary  to 
cluster  grid  points  In  the  leading  and  trailing 
edge  regions.  Minimum  values  of  As  «  0.002  at 
the  leading  edge  and  0.005  at  the  trailing  edge 
were  necessary  to  resolve  the  riow  gradients  near 
the  nose  and  to  attain  regular  behavior  at  the 
trailing  edge.  The  As  spacing  was  Increased  to 
approximately  0.025  at  midchord.  This  spacing 
caused  the  computed  shock  waves  to  be  smeared 
somewhat,  but  the  aerodynamic  coefficients  were 
affected  little  by  this. 


Numerical  Algorithm  and  Characteristics 

The  numerical  method  used  to  solve  Eq.  !’)  is 
an  alternating-direction  fully  implicit  :ADi;. 
approximate-factorization  scheme,  described  In 
detail  in  Refs.  5-7.  We  used  the  Method  3  option, 
l.e.,  scalar  pentadlagonal  inversion,^  of  ARC2D- 
Version  1.50.  The  finite-difference  approxima¬ 
tions  to  Eq.  (1)  are  solved  as  a  series  of  time¬ 
like  iterations,  with  the  steady  state  approached 
asymptotically.  However,  for  computational  effi¬ 
ciency,  the  time  step  At  is  varied  throughout 
the  computational  domain  and  with  iteration 
number,  so  that  the  Intermediate  results  are  not 
"time  accurate."  Thus  this  variable  time  step  is 
a  numerical  parameter  that  can  be  varied  by  the 
user  of  the  code  to  enhance  numerical  stability 
and  to  accelerate  the  convergence. 

Another  numerical  parameter  that  can  be 
varied  by  the  user  Is  the  artificial,  or  numeri¬ 
cal,  dissipation.  Fourth-order  implicit  and 
explicit  nonlinear  dissipation  terms  are  added  to 
Improve  the  nonlinear  stability  of  the  algorithm 
and  to  enhance  Its  convergence,  especially  for 
transonic  flows  with  shock  waves. 

We  found  by  numerical  experimentation  that 
values  of  At  and  the  dissipation  coefficients 
close  to  the  default  values  were  satisfactory  in 
most  cases,  although  the  oases  with  the  largest 
separated  zones  required  special  treatment.  Also, 
the  grid  spacing  at  the  body  was  altered  for  the 
cases  at  the  extreme  values  of  Reynolds  number. 
However,  the  convergence  histories  shown  in  Fig.  k 
are  representative  of  most  of  the  results  given  in 
this  paper.  The  engineering  quantities  of  Inter¬ 
est  (l.e.,  lift,  drag,  and  pitching-moment  coeffi¬ 
cients)  normally  converged  to  within  IS  of  their 
final  values  within  about  1200  Iterations,  or 
approximately  5  min  CPU  time  on  the  Ames  Cray  XM-P 
supercomputer .  However,  we  always  carried  the 
calculations  well  beyond  that  point,  until  the 
residuals  had  dropped  another  order  of  magnitude, 
as  Indicated  in  Fig.  k.  The  CPU  times  were  about 
0.31  sec/lteratlon  for  the  193  *  6k  grid. 

Finally,  we  have  developed  some  subjective 
Impressions  of  the  overall  accuracy  of  the  com¬ 
puted  results,  which  are  based  on  our  extensive 
numerical  experiments  and  on  comparisons  with  wind 
tunnel  data,  as  discussed  in  the  following  sec¬ 
tion.  As  noted  above,  the  calculations  normally 
converged  to  engineering  accuracy,  and  also  to 
within  the  limits  of  the  numerical  errors,  in 
about  1200  iterations.  Apart  from  cases  with  mas¬ 
sive  separation,  the  results  for  lift  appear  to  be 
less  than  5t  in  error;  as  shown  below,  the  calcu¬ 
lations  match  or  slightly  overpredict  the  measured 
values  of  CL.  The  absolute  uncertainty  in  the 
calculated  drag  coefficient  Is  higher,  because  of 
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airfoil  from  many  sources  nave  seen  screened  ana 
correlated,  and  those  data  from  Refs.  ’  and  "3-2 
have  been  assembled  in  Fig.  5.  These  sources 
include  data  that  are  considered  to  be  virtually 
interference-free, results  which  have  been 
polished  with  the  latest  state  of  the  art  ir.  wal 
correction  techniques,17  a  recent  set  of  details 
measurements  of  exceptional  quality,  but  wnic.n 
require  significant  angle-of -attack  corrections, 
and  a  benchmark  data  set  from  wartime  Germany.  i 
Also  shown  are  representative  results  that  are 
commonly  used  in  the  helicopter  indus- 
try . 1 ■ 1 5,20-22  The  Reynolds  numbers  vary 


•  Gregory  and  Wilby,  NPL  36  in.  •  m  in. 
trip 

9  -  Vidal  et  al.,1®  Calspan  8  ft;  trip 


9  •  Harris,10  Langley  8  ft  TPT;  trip 


1  800  1000  1800  2000  2800  i 
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lJio-’o 

3000 


□  •  Noonan  and  Bingham,  0  Langley 
6  in.  «  28  in. ;  no  trip 

9  •  Noonan  and  Bingham,20  Langley 
6  in.  ■  28  in.;  trip 

A  -  Prouty,1  Lockheed  15  in,  «  A8  in.  2-D 


Fig.  M  Convergence  history  of  a  typical 
calculation. 


-  Cumber t  and  Newman.1'  Langley  0.3~m 
Cryogenic 


the  sensitivity  to  the  grid,  to  the  turbulence 
model,  and  to  the  transition  location.  The  latter 
could  account  for  up  to  ±0.0020  in  CD  below  drag 
divergence;  otherwise,  the  errors  are  probably  on 
the  order  of  ±10%  over  the  range  of  conditions 
considered  herein.  Our  estimates  of  the  uncer¬ 
tainty  in  the  pitching  moment  coefficient  are  even 
more  speculative:  the  greater  of  either 
Cm  -  ±0.005  or  ±101.  While  these  uncertainty 
estimates  are  comparable  to  those  of  the  best 
modern  wind  tunnels,  they  are  not  as  low  as  the 
aircraft  Industry  would  like.11* 


Results  and  Discussion 


As  indicated  in  Fig.  2,  approximately  90 
combinations  of  airfoil,  Mach  number  and  angle  cf 
attack  were  computed,  mostly  at  Re  •  6  «  10  . 

The  results  are  summarized  in  this  section. 


7  -  Gothert,1*  DVL  2.7-m 
O  •  Jepson,21  Tanner,22  Llzak,15,  UTRC  8  ft 
©  »  Jepson,21  DTNSRDC  7  ft  «  10  ft 
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Experimental  Data  for  Comparison  and  Validation 


As  noted  earlier  in  connection  with  Table  1, 
wind  tunnel  results  differ,  according  to  the  par¬ 
ticular  characteristics  of  each  facility.  For  the 
purposes  of  this  paper,  data  for  the  NACA  0012 


2  2  .4  &  .8  .7  .8  .9  1.0  1.1 

MACH  NUMBER.  M„ 

Fig.  5  Experimental  results  for  the  NACA  0012 
airfoil. 


approximately  ’0®  to  10^,  and  transition  is  Both 
'ixed  and  free. 

In  fig.  5  and  in  several  subsequent  figures, 
the  shaded  Bands  are  our  assessments  and  interpre¬ 
tations  of  the  limits  of  this  aggregate  of  data; 
that  is,  the  extent  to  which  one  may  reasonably 
claim  to  know  the  force  and  moment  characteristics 
of  the  NACA  0012  profile  in  practice.  Seme  data 
have  been  discounted,  for  example,  the  low  values 
of  from  Refs.  1  and  21,  which  are  probably 

contaminated  by  wall-interference  effects.  Also, 
the  Calspan16  values  of  drag  coefficient  at  low 
Mach  number  appear  to  be  too  high,  perhaps  due  to 
the  type  of  boundary  layer  trip  that  was  used.  In 
the  rest  of  this  paper,  our  calculated  results 
will  be  compared  with  experimental  data  bands  that 
have  been  screened  from  Refs.  1  and  15-22  in  this 
manner.  It  must  be  emphasized  that  the  validity 
of  the  calculations  can  only  be  established  to 
tolerances  within  these  limits,  and  that  no  single 
experiment  exists  that  Is  satisfactory  for  code 
validation. 

Results  for  the  NACA  0012 

The  numerical  results  for  the  NACA  0012  pro¬ 
file  as  a  function  of  Mach  number  and  angle  of 
attack  are  shown  in  Figs.  6-10.  As  noted  In  the 
description  of  the  numerical  method,  these  calcu¬ 
lations  were  all  performed  at  Re  -  6  «  10®,  and 
boundary  layer  transition  was  assumed  to  occur 
at  x/c  -  0.01 .  Figures  7-9  Include  the  experi¬ 
mental  data  described  above;  also  shown  for  refer¬ 
ence  are  the  results  of  a  standard  transonic 
small-disturbance  lnvlscld  code,  LTRAN2. 

Strong  nonlinear  compressibility  effects  are 
evident  In  CL  and  CQ  as  Ha  increases;  e.g., 
the  well-known  drag  divergence  at  Ma  -  0.78  and 
the  abrupt  reduction  In  lift  for  Ma  >  0.82.  The 
spread  in  (L/D)m#Jt  is  much  greater,  but  a 
substantial  reduction  with  Increasing  Ma  is 
apparent.  The  present  results  tend  to  fall  below 
the  middle  of  the  experimental  band,  probably 
because  the  transition  point  was  chosen  at 
x/c  •  0.01 . 

It  is  clear  that  the  airloads  computed  by  the 
present  ARC2D  viscous  code  exhibit  all  the  fea¬ 
tures  that  have  been  revealed  by  the  experiments, 
and  the  values  of  the  coefficients  are  approxi¬ 
mately  within  the  ranges  of  the  measured  ones. 
There  are  two  relatively  minor  discrepancies: 
first,  the  computed  values  of  at  subsonic 
Mach  numbers  tend  to  be  slightly  higher  than  the 
data;  and  second,  the  computed  drag  coefficients 
follow  the  upper  extremity  of  the  experimental 
band.  On  the  other  hand,  the  nonlinear  lnvlscld 
computations  predict  even  higher  at  subsonic 
Mach  numbers,  and  the  lnvlscld. results  are 


0  123486789  10 


ANGLE  OF  ATTACK,  a 

Fig.  6  Calculated  results  for  the  NACA  0012 
airfoil  with  trip,  a)  Lift  coefficient  vs  angle 
of  attack;  b)  drag  coefficient  vs  angle  or  attack. 

completely  unsatisfactory  in  the  range 
0.8  <  Ma  <  0.9.  This  Is  associated  with  the  shock 
wave-boundary  layer  interaction  that  dominates  the 
real  flow,  as  shown  below. 

As  is  well  known,  the  drag  divergence  is 
related  to  the  formation  of  shock  waves  on  the 
airfoil  for  Ma  >  0.78.  Pressure  distributions 
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O  TRANSONIC  SMALL  DISTURBANCE 
•  NAVIER-STOKES 
TP77X  EXPERIMENT 


MACH  NUMBER. 


Fig.  8  Comparison  of  drag  at  zero  lift  for  the 
NACA  0012  airfoil. 


and  boundary-layer  Information  that  help  explain 
the  highly  nonlinear  behavior  of  CL  in  Fig.  7 
are  shown  in  Fig.  10,  corresponding  to  Points  A-D 
in  Fig.  7.  Figure  10  illustrates  the  buildup  of 
shock  strength  and  shock  wave/boundary  layer 
interaction  with  increasing  Mach  number. 

First,  a  small  separation  bubble  appears  cn 
the  upper  surface  at  Mw  -  0.80,  and  although  the 
lift  is  greater  than  predicted  by  linear  theory, 
it  is  significantly  below  the  nonlinear  lnvlscld 
value.  At  ■  0.85,  the  interaction  between  the 
shock  wave  and  the  boundary  layer  produces  exten¬ 
sive  separation  on  the  upper  surface  and  a  small 
separation  bubble  on  the  lower  surface,  causing  a 
further  reduction  in  lift-curve  slope.  This  trend 
continues  up  to  •  0.88  where,  according  to 
the  viscous  calculations,  a  minimum  value  of 
dC^/da  ■  0  occurs.  At  this  point,  the  upper- 
surface  shock  wave  is  stronger  than  the  lower,  but 
the  separated  viscous  layer  pushes  it  farther 
upstream;  consequently,  the  net  lift  is  approxi¬ 
mately  zero,  but  the  pitching  moment  derivative 
(not  shown)  attains  a  maximum  positive  value  of 
dC^/da  *  0.04/deg. 


At  higher  Mach  numbers,  the  shock  waves  on 
both  surfaces  move  to  the  trailing  edge,  the  vis¬ 
cous  effects  become  relatively  unimportant,  and 
the  rates  of  change  of  CL,  CD,  and  CM  with  mb 
become  much  less.  No  further  significant  qualita¬ 
tive  changes  occur  up  to  the  maximum  Mach  number 
considered,  MB  ■  1.10.  It  is  interesting  to  note 
that  the  lnvlscld  code  becomes  surprisingly  accu¬ 
rate  again  for  M  >  0.9<t. 


.8 


1.2 


X/C 


Fig.  10  Pressure  distributions  for  the  NACA  0012  slrfoll;  a  ■  1/2°,  Re  -  6  *  10^,  S  denotes  boundary- 

layer  separation,  R  reattachment. 


Correlation  of  Transonic  Drag  Characteristics 


Recently  F.  D.  Harris  (private  communication) 
reminded  the  authors  of  Sprelter's2^  work  with 
similarity  parameters  as  an  aid  In  correlating  the 
drag  due  to  shock  waves  In  the  transonic  regime. 
The  basic  transonic  similarity  parameters  are 


c.  m^/3(y  ♦  1)1/3 


correlate  approximately  uniquely  with  M  and  the 
position  of  the  maximum  airfoil  thickness. 

Figure  1 1  shows  that  our  calculated  results 
of  CD  versus  M  correlate  well  for  the  four 
airfoils  considered.  Also  shown  In  the  figure  are 
Harris'  correlation  of  a  number  of  finite  wing 
tests,  and  the  band  of  data  from  Refs.  1  and 
15-22.  As  noted  earlier  In  connection  with 
Fig.  8,  our  numerical  results  follow  the  upper 
limit  of  the  available  measurements. 


M2  -  1 


My3(t/c)2/3(Y  *  ni/3 


where  C0^  l9  increment  In  drag  coefficient 

above  the  subsonic  level,  y  is  the  specific  heat 
ratio,  and  t/c  Is  the  airfoil  thickness  ratio. 
Harris  (private  communication)  has  found  to 


Results  for  the  Vertol  VR-8  Airfoil 

The  transonic  pitching  moment  behavior  has 
particularly  Important  consequences  on  the 
advancing  blade  of  a  rotor  In  high-speed  forward 
flight.  Large  torsional  loads  can  be  created  by 
the  rapid  growth  in  that  occurs  when  the 
aerodynamic  center  of  pressure  shifts  rearward  In 
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Fig.  11  Correlation  of  drag  results  In  transonic 
similarity  parameters. 

the  transonic  regime,  a  phenomenon  known  as  "Mach 
tuck." 

Figures  12  and  13  show  the  computed  pitching 
moment  behavior  and  drag  coefficient  of  the  VR-8 
section  at  zero  lift  compared  with  wind  tunnel 
measurements.1^  The  calculations  were  performed 
both  at  a  fixed  Reynolds  number,  corresponding  to 
the  NACA  symmetrical -airfoil  results  described 
above,  and  at  the  varying  Reynolds  numbers  of  the 
experiment,  8-13  «  10°.  Transition  was  assumed 
at  x/c  -0.01  in  all  cases.  As  shown  in  the 
figures,  the  results  are  insensitive  to  Reynolds 
number  In  this  limited  range. 

The  large  Mach-tuck  dip  that  occurs  for 
0.8  <  <  0.9  is  clearly  evident  in  Fig.  12. 

The  calculated  results  shown  In  Figs.  12  and  13 
are  in  excellent  agreement  with  the  measurements, 
except  for  the  maximum  negative  value  or  and 
the  Mach  number  at  which  it  occurs.  The  values  of 
(L/D)max  at  Ma  •  0.60  and  CQ  versus  M  were 
snown  earlier  in  Figs.  9  and  11,  respectively.  In 
both  cases,  the  aerodynamic  performance  Is  supe¬ 
rior  to  that  of  the  NACA  0012. 


- -O  NAVIER-STOKES  R.  -  6  <  106 

□  NAVIER-STOKES.  EXPERIMENTAL  R. 


Fig.  12  Pitching  moment  vs  Mach  number  for  the 
VR-8  airfoil  at  CL  -  0. 


Fig.  13  Drag  coefficient  vs  Mach  number  for  the 
VR-8  airfoil  at  •  0. 

and  shock  wave/boundary  layer  Interaction  with 
Increasing  Mach  number.  However,  the  actual  pres¬ 
sure  distributions  are  considerably  different  from 
those  shown  In  Fig.  10.  The  cases  between 
Mb  •  0.85  and  0.90,  which  encompass  the  peak  in 
-CM,  coincide  with  the  maximum  amount  of  shock- 
induced  separation  on  the  airfoil.  The  disagree¬ 
ment  with  the  experimental  data  for 
0.85  <  Mw  <  0.90  could  be  due  to  errors  in 


The  pressure  distributions  and  boundary-layer 
Information  relevant  to  the  nonlinear  behavior 
shown  In  Figs.  12  and  13  are  given  in  Fig.  10.  As 
In  the  NACA  0012  results  considered  earlier,  the 
figure  Illustrates  the  buildup  of  shock  strength 


modeling  the  boundary  layer  transition,  as  the 
experimental  transition  characteristics  were  not 
specified.  On  the  other  hand,  the  spread  in  the 
available  data  on  the  NACA  0012  airfoil  (e.g.. 
Figs.  5,  7,  8,  and  9)  suggests  that  it  is  unwise 
to  place  too  much  faith  in  a  single  experiment. 


Fig.  1U  Pressure  distributions  for  the  VR-8  airfoil  at  CL  •  0;  S  denotes  boundary-layer  separation, 

R  reattachment. 


Effects  of  Reynolds  Number 


Subject  to  the  limitations  of  the  empirical 
models  for  turbulence  and  transition,  the  computa¬ 
tions  can  be  run  over  virtually  an  unlimited  range 
of  Reynolds  numbers.  Representative  subsonic  and 
transonic  results  for  the  NACA  0012  airfoil 
between  1  and  100  *  10°  are  shown  In  Fig.  15  for 
transition  fixed  at  x/c  ■  0.01.  The  results  show 
a  monotonic  decrease  In  CD  with  Increasing  Re, 
whereas  the  lift  monotonically  Increases. 

The  pressure  distributions  corresponding  to 
the  Ms  •  0.80  case  are  given  In  Fig.  16.  These 
results  are  all  qualitatively  the  same,  but  the 
shocx-induced  boundary-layer  separation  Is  a 
strong  function  of  Reynolds  number.  For  example, 
tne  separation  zone  on  the  upper  surface  at 
Re  *  10°  extends  from  the  shock  wave  to  the 
trailing  edge,  whereas  It  Is  only  a  short  bubble 
at  Re  •  10',  and  there  is  no  separation  at 
Re  •  10s.  It  Is  interesting  to  note  that  although 
the  subsonic  drag  coefficients  are  considerably 


REYNOLDS  NUMBER,  Re 


Fig.  15  Effect  of  Reynolds  number  on  lift  and 
drags  NACA  0012,  M.  •  0.5  and  0.8,  a  -  i*. 
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Fig.  16  Pressure  distributions  on  the  NACA  0012  airfoil  at  M  ■  0.80  and  a  -  1*. 


less  than  the  transonic  values,  the  percentage 
change  with  Re  is  actually  greater. 


SugmaKy  and  Conclusions 

In  this  paper,  we  have  applied  a  new  aerody¬ 
namic  tool  to  the  study  of  helicopter  airfoil 
characteristics.  We  have  shown  that  the  computed 
airloads  reproduce  completely  the  experimental 
behavior  of  representative  airfoils  across  the 
transonic  regime.  In  addition,  the  computational 
details  of  the  flow  fields,  the  surface  pressure 
distributions,  and  the  viscous-layer  characteris¬ 
tics  enable  us  to  trace  the  evolution  of  the  phys¬ 
ical  changes  tnat  occur  as  Ma  or  Re  increases. 
Descriptions  of  the  complicated  development  of 
shock  waves,  shock  wave-boundary  layer  interac¬ 
tion,  and  shock-induced  separation  supplement  the 
information  that  has  been  obtained  heretofore  in 
wind  tunnels. 

In  validating  our  calculations  and  assessing 
the  accuracy  of  the  results,  including  extensive 


grid-refinement  studies  and  comparisons  with  data 
from  numerous  wind  tunnels,  we  have  defined  more 
precisely  the  capabilities  and  limitations  of  the 
code  ARC2D.  This  important  aspect  of  the  investi¬ 
gation  has  also  revealed  ways  that  the  computa¬ 
tions  can  complement  wind-tunnel  tests,  by  provid¬ 
ing  flow-field  details  that  are  difficult  to 
measure  and  by  extending  the  range  of  flow  param¬ 
eters  beyond  the  capabilities  of  existing  wind 
tunnels.  The  code  has  now  progressed  from  a 
purely  research  stage  to  almost  a  production 
stage,  where  it  can  be  run  by  specialists  in  the 
helicopter  industry, 

He  claim  that  the  numerical  results  presented 
in  this  paper  are  as  accurate  as  contemporary 
experimental  measurements,  and  that  new  physical 
insights  into  transonic  airfoil  behavior  can  be 
obtained  from  the  details  of  these  calculations. 
However,  we  emphasize  that  these  and  similar 
emerging  CFD  capabilities  represent  new  tools  to 
complement  experimental  measurements,  rather  than 
substitutes  for  wind  tunnel  testing.  Each 
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approach  has  Its  unique  capabilities,  limitations , 
economics,  and  manpower  requirements.  The  future 
of  helicopter  technology  Is  best  served  by  an 
intelligent  combination  of  the  two. 
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